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Abstract. We review the collider phenomenology of neutrino physics and the 
synergetic aspects at energy, intensity and cosmic frontiers to test the new physics 
behind the neutrino mass mechanism. In particular, we focus on seesaw models 
within the minimal setup as well as with extended gauge and/or Higgs sectors, and 
on supersymmetric neutrino mass models with seesaw mechanism and with i?-parity 
violation. In the simplest Type-I seesaw scenario with sterile neutrinos, we summarize 
and update the current experimental constraints on the sterile neutrino mass and its 
mixing with the active neutrinos. We also discuss the future experimental prospects 
of testing the seesaw mechanism at colliders and in related low-energy searches for 
rare processes, such as lepton flavor violation and neutrinoless double beta decay. The 
implications of the discovery of lepton number violation at the LHC for leptogenesis 
are also studied. 


1. Introduction 

With the discovery of a Higgs boson at the Large Hadron Collider (LHC) 
having properties consistent with the Standard Model (SM) expectations [^, we are 
tantalizingly close to verifying the Higgs mechanism responsible for the SM gauge boson 
and charged fermion masses. What remains missing though is an understanding of the 
neutrino masses. The observation of neutrino oscillations in solar, atmospheric, reactor 
and accelerator neutrino data demonstrates that at least two of the three active 
neutrinos have a non-zero mass and that individual lepton flavor is violated. This 
provides the hrst and so far only conclusive experimental evidence for the existence of 
New Physics beyond the SM. 

In the SM, neutrinos have only one helicity state up, and therefore, cannot acquire 
a Dirac mass, unlike the charged fermions, after the electroweak symmetry breaking 
(EWSB) by the vacuum expectation value (VEV) of the SM Higgs held. Just adding 
‘by hand’ a right-handed (RH) neutrino held N per generation to the SM, one could in 
principle generate a Dirac mass term for neutrinos. However, to get sub-eV left-handed 
(LH) neutrino masses as required by the neutrino oscillation data, one needs the Dirac 
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Yukawa couplings to be extremely tiny < 10“^^. There is no theoretical justihcation 
for the large disparity between such small neutrino Yukawa couplings and other SM 
Yukawa couplings. Moreover, such a scenario would be rather uninteresting from an 
experimental point of view. Therefore, we will take the more optimistic viewpoint that 
some other New Physics might be responsible for the observed smallness of neutrino 
masses. 

Being electrically neutral, neutrinos can in principle have a Majorana mass term 
of the form where C being the charge conjugation matrix, and 

I = e, /i, r is the lepton flavor index. However, since pu is part of the SU (2)^ 
doublet field with lepton number L = +1, the above Majorana mass 

term transforms as an SU{2)l triplet, i.e. it is not gauge invariant, apart from breaking 
the global L and B — L symmetries of the SM by two units. On the other hand, the only 
known source of lepton number violation (LNV) in the SM is via non-perturbative weak 
instanton effects through the Adler-Bell-Jackiw anomaly |^|^. These non-perturbative 
effects only violate the B + L symmetry |^, but preserve the B — L symmetry to all 
orders. Therefore, neutrino masses cannot be induced even by non-perturbative effects 
in the SM. 

One might wonder whether quantum gravity effects might be sufficient to explain 
the tiny neutrino masses. As long as gravity is treated perturbatively and respects the 
B — L symmetry of the SM, neutrinos remain massless. However, non-perturbative 
gravitational effects, e.g. black holes and worm holes, do not respect global symmetries 
and could induce non-zero neutrino masses in the effective low energy Lagrangian. In 
the context of the SM, these effects are at most of order n^/Mpi 10 ® eV (where 
V is the SM Higgs VEV and Mpi is the Planck mass), which is still three orders of 
magnitude below that required to satisfy the atmospheric neutrino data [^. 

Therefore, it seems more natural to invoke some {B — L)-violating New Physics 
beyond the SM at a scale A <C Mpi to explain the observed neutrino masses (^. 
Within the SM, these {B — L)-breaking effects can be parametrized through an effective 
dimension-5 Weinberg operator of the form Xu/LgLf^^/A |^, where $ = {(p'^, is 
the SM Higgs doublet. It is easy to see that there are only three ways to obtain this 
Weinberg operator at tree-level, using only renormalizable interactions, i.e.. 


The product of L^ and <I> forms a fermion singlet: {LJ^){LJ,^)/A. This is widely 
known as the Type-1 seesaw mechanism ll-Tsl. Here the intermediate heavy 


particles are clearly fermion singlets, identified as the RH Majorana neutrinos Njia- 
The product of Li and L^/ forms a scalar triplet: {LjaaLir){^'^aa^)/A, where Ua’s 


111 


are the usual Pauli matrices. This is known as the Type-II seesaw mechanism |15f - 
19 . Here the intermediate heavy particle is a scalar triplet A = (A’*"’*', A+, A°). 

The product of Li and <I> forms a fermion triplet: {Lja°'^){Lj,a'^^)/A. This is 
known as the Type-Ill seesaw mechanism 20 . Here the intermediate heavy particle 
is a fermion triplet S = (E'*', E“). 


One could also construct a hybrid seesaw model using more than one of these different 
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types of seesaw mechanism. Non-minimal variations of the seesaw mechanism with 


higher multiplets have also been constructed; see e.g. 21,22 


The origin of the bare Majorana mass terms responsible for the explicit B — L 
violation can be understood from natural implementations of the seesaw mechanism 
in ultra-violet (UV) complete theories, e.g. in the Left-Right Symmetric Model 
(LRSM) |2§g and in SU{b) |^, ^O(IO) and Eg |29}|^ Grand Unihed 

Theories (GUTs); for a review, see e.g. (3^ . 

For the conventional seesaw scenarios mentioned above, the light neutrino masses 
are inversely proportional to a large lepton-number breaking scale (hence the name 
‘seesaw’). There exists an important variation, where the light neutrino masses are 
directly proportional to a small lepton-number breaking scale. This is known as the 
inverse seesaw mechanism (see also 29 ,^) and can be regarded as an extension of 


the Type-I seesaw, since we introduce two sets of SM singlet fermions, instead of one. 
Experimentally, the main distinguishing feature of this scenario is the pseudo-Dirac 
nature of the heavy SM-singlet fermions, in contrast with the purely Majorana nature 
of the singlet fermions in the Type-I seesaw scenario. 

Another interesting class of neutrino mass models uses the radiative mechanism 
which, unlike the tree-level seesaw models, can go beyond the effective dimension-5 
operator and generate tiny neutrino masses at an arbitrary loop-level. Some of the 
simplest and predictive radiative seesaw models are given in 33 43 . In supersymmetric 


theories, there exists yet another possibility of generating non-zero neutrino masses 
without the need of a seesaw mechanism. The minimal way to incorporate {B — L)- 
breaking with the MSSM particle content is by bilinear i?-parity violation |44|-[5^. For 


a review on various low-scale neutrino mass models, see e.g. 54 


In a bottom-up approach, the scale of new physics responsible for {B — L)-breaking 
is a priori unknown. Since this review is about the collider tests of neutrino physics, we 
will only consider those scenarios with a {B — L)-breaking scale < (9(TeV) accessible 
to foreseeable collider experiments. The rest of the review is organized as follows: in 
Section we discuss the collider prospects of sterile neutrinos appearing in the minimal 
Type-I seesaw and its singlet extensions within the SM gauge group. In Section 
we discuss the collider prospects of minimal Type-II and Type-III seesaw models. In 
Section]^ we review some seesaw models with extended gauge sectors, namely, with an 
additional U{1) in Section 4.1| and with the Left-Right (L-R) Symmetric gauge group 
in Section 4.2[ In Section we discuss the signihcance of the observation of LNV for 
leptogenesis. In Section we briefly discuss the supersymmetric neutrino mass models. 
Finally, in Section we present our conclusions and future prospects. 
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The simplest renormalizable extension of the SM for nnderstanding the smallness of the 
LH nentrino masses is dehned by the interaction Lagrangian 


-Cy — hiaLi^Njia + H.C. , 


( 1 ) 


where = ia 2 ^* and Nji^ are SM singlet nentral fermions, also known as the sterile 
nentrinos, since they cannot directly participate in the SM charged-cnrrent (CC) and 
nentral-cnrrent (NC) interactions in the absence of any mixing with the active nentrino 
sector. In ([^, a = 1,2, ■■■ ,J\f is the sterile nentrino flavor index and are the 
dimensionless complex Yukawa couplings. From the structure of Q, we see that 
the fermions must necessarily be right-chiral; hence, they are also known as RH 


neutrinos. This distinction will be naturally justihed in the LRSM (see Section 4.2) 


Assuming that ([^ is the only source of neutrino masses and oscillations, we need at 
least two or three RH neutrinos, depending on whether the lightest active neutrino is 
massless or not||] 

In the Higgs phase after EWSB, the term ([^ generates a Dirac mass Md = hv. 
Since the RH neutrinos carry no SM gauge charges, one can also write a Majorana mass 
term 


— = -^iMN)ayNRa^Ry + 


( 2 ) 


while preserving gauge invariance. The term ([^ implies that the hypercharge of Nr^ is 
zero, and therefore, from ([^, we deduce that the hypercharges of the lepton and Higgs 
doublets are the same. Thus, the requirement of cancellation of gauge chiral anomalies 
implies charge quantization, provided that the neutrino mass eigenstates are Majorana 
helds |5§]^|| 

The terms ([^ and (|^ together lead to the following neutrino mass matrix in the 
flavor basis ^Ra}- 


M,. = 


0 Md 

MJ) Mj\f 


( 3 ) 


For 1 (with ||M|| = being the norm of matrix M), the light 

neutrino masses and mixing are given by the diagonalization of the effective mass matrix 


and the left-right neutrino mixing parameter is given by 
Type-I seesaw mechanism 11-15 , as mentioned in Section]^ 


( 4 ) 


MrM^^. This is the 


f The current upper limits on the absolute active neutrino mass scale from the kinematics of tritium 


/3-decay 55,56 , as well as the cosmological bounds on the sum of neutrino masses 57 still allow for 


three non-zero active neutrinos. 

§ This is true regardless of the number of generations, in contrast with the SM (without RH neutrinos) 


where anomaly cancellation implies charge quantization only in the one generation case 62 63 . 
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From the above discussion, it is clear that there are two key aspects of the Type-I 
seesaw mechanism that can be probed experimentally: the Majorana mass M^v of the 
mostly sterile neutrinos and their mixing Vin with the active neutrinos. The Majorana 
nature of both light and heavy neutrinos can in principle be tested via the classic 


LNV process of neutrinoless double beta decay for reviews, see e.g. 64,65 


However, an observation of 0z//9/9 does not necessarily provide us with information on 
the mixing whose effects on the 0z//?/9 amplitude may be sub-dominant, compared 
to purely left- or right-chiral contributions. On the other hand, if the mixing effects are 
non-negligible, they could be inferred from non-unitarity of the light neutrino mixing 
matrix 6^68 , as well as in low-energy observables for lepton flavor violation (LFV), 
lepton non-universality and electroweak precision tests 68-73 . However, these low- 


energy observables by themselves do not prove the Majorana nature of heavy neutrinos, 
since models with pseudo-Dirac heavy neutrinos can also yield large non-unitarity and 
LFV effects 74-78 . As we will discuss in this section, the collider experiments at 


the energy frontier provide a simultaneous probe of both the aspects of seesaw, if the 
heavy neutrinos are kinematically accessible. This is complementary to the low-energy 
searches of LNV and LFV at the intensity frontier. 

In a bottom-up approach, the RH neutrinos can just be introduced ‘by hand’ as the 
only new particles beyond the SM, as e.g. in the Neutrino SM (i/SM) [TO , 80 . In this 


case, the RH neutrino masses in ([^ are largely unconstrained, even with the requirement 
of satisfying the neutrino oscillation data. For instance, if the Yukawa couplings 
\hia\ ~ one expects from Q that the Majorana mass scale Mjv ~ 0(10^^ GeV), 

which is close to the upper limit due to partial-wave unitarity: A < 47rn^/\/3My [^ . 
For any smaller value of Mj^, one could hnd the associated values of Yukawa couplings, 
as required to £t the light neutrino data using the seesaw formula (j^. In fact, there exist 
seesaw models with spanning over a wide range of scales, i.e. from eV to GUT scale 
(for a review, see [^). However, from the experimental point of view, only the scenarios 
with Misf < C>(TeV) offer the possibility of being tested in foreseeable future. Even from 
a theoretical point of view, the naturalness requirement that electroweak corrections to 
the bilinear Higgs mass operator should not exceed 6fi ~ 0{1 TeV) imposes 

an upper bound of Mjv 4 x 10^ GeV on the lightest RH neutrino mass |^-86 


Supergravity models of inflation impose an additional upper bound on the reheating 
temperature Tr < 10® — 10® GeV (and hence, on the RH neutrino masses), as required 
to avoid overproduction of gravitinos whose late decays may otherwise spoil the success 
of Big Bang Nucleosynthesis (BBN) [^ 93 . Another reason for considering the low- 


scale seesaw mechanism is that it provides a unique opportunity to test the link between 
the neutrino mass mechanism and the observed matter-antimatter asymmetry in our 
Universe via the mechanism of leptogenesis (for a review, see e.g. [^). Motivated by 
all these theoretical and experimental considerations, we will focus on the low-scale 
seesaw scenarios. 

We should also mention here that there has been a lot of recent interest in GeV- 


scale seesaw models, especially within the i/SM framework 79,80 , which can explain not 
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only neutrino masses, but also matter-antimatter asymmetry and Dark Matter (DM), 
in a minimal setup (for a review, see e.g. [^). This scenario relies on the possibility 
that two sterile neutrinos with masses typically in the range of 1-10 GeV, well below 
the critical temperature of the electroweak phase transition, lead to an enhanced lepton 


asymmetry via coherent CP-asymmetric oscillations 80,96-99 , which is of 0{h^) but 
could be sufficient to explain the observed matter-antimatter asymmetry. This requires 
a very specihc choice of model parameters, but might be accessible to current and 
future experiments at the intensity and cosmic frontiers. In particular, the lightest RH 
neutrino could play the role of a DM, if its mass is in the keV range (for a review, see 
e.g. [Too]). This has gained further attention in light of the recent observations of an 


unidentified X-ray line at energy of 3.5 keV 101-103 , which could be explained as due 


to the radiative decay of a sterile neutrino DM with mass of 7 keV 
In an ultra-violet complete theory, such as in the LRSM 
GUT 


23-25 and in S'O(IO) 


104 , the RH neutrinos constitute an integral part of the particle spectrum, as 


required by anomaly cancellation, and the RH neutrino mass scale is no longer an adhoc 
parameter, but intimately related to the P(l)B-L-breaking scale. In the S'O(10)-type 
theories, quark-lepton unihcation implies that the Yukawa couplings appearing in (j^ 
are of similar order of magnitude as the up-quark Yukawa couplings, which means that 
the seesaw scale Mjy ~ 10^“^ GeV. Thus, even though the GUT seesaw models are quite 


elegant and predictive 28 , testing their ultra-heavy particle spectrum experimentally, 
going beyond the ‘grand’ desert, is a formidable task. On the other hand, simpler 
gauge-extended seesaw models, such as the L-R seesaw, could easily accommodate 
an experimentally-accessible RH neutrino mass scale. We will discuss some of these 
Section ^ 


scenarios m 


2.1. Low-Scale Singlet Seesaw Models 


In the traditional “vanilla” seesaw mechanism 11 15 , the left-right neutrino mixing is 
given by 


VfM ~ 




M, 


< 

rs_/ 


10 


-6 


N 


100 GeV 


M, 


N 


due to the smallness of the light neutrino mass < 0.1 eV 57 . 


( 5 ) 


Thus for a low 


seesaw scale in the sub-TeV to TeV range, the experimental effects of the light-heavy 
neutrino mixing are expected to be too small, unless the RH neutrinos have additional 
interactions, e.g. when they are charged under U{1)b-l- However, there exists a class 


of minimal SM plus low-scale Type-I seesaw scenarios 36,119-130 , where Vin can be 


sizable while still satisfying the light neutrino data. This is made possible by assigning 
II An alternative low-scale LNV can be realized in models with spontaneous global B — L violation, 


which implies the existence of a massless Goldstone boson called the Majoron 105 


interesting phenomenological consequences for Ovjdfi 108 -112 , leptogenesis 113 


107 with potentially 


DM 114-116 and 


inflation 117 . In the context of supersymmetry, both B — L and it-parity can be spontaneously broken 


by the VEV of the RH sneutrino field, leading to a minimal model of i?-parity violation 118 . 
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specific textures to the Dirac and Majorana mass matrices in the seesaw formula Q. The 
stability of these textures can in principle be guaranteed by enforcing some symmetries 
in the lepton sector 121, 122, 128,131,1132 . We will generically assume this to be 


the case for our subsequent discussion on the collider signatures of low-scale minimal 
seesaw, without referring to any particular texture or model-building aspects. Also, 
unless otherwise specihed, we will use a model-independent phenomenological approach, 
parametrized by a single heavy neutrino mass scale Mjv and a single flavor light-heavy 
neutrino mixing Wjv, assuming that the mixing effects in other flavors P ^ I are sub¬ 
dominant. Although this assumption may not be strictly valid for a realistic seesaw 
model satisfying the observed neutrino oscillation data, it enables us to derive generic 
bounds on the mixing parameter, which could be translated or scaled appropriately in 
the context of particular neutrino mass models (see e.g. 133| ). 

Another natural realization of a low-scale seesaw scenario with large light-heavy 


neutrino mixing is the inverse seesaw model 31 , where one introduces two sets of SM 


singlet fermions {Npia, Sip} with opposite lepton numbers, i.e. T(A'h) = -|-1 = —L{Sl) 
In this case, the neutrino Yukawa sector of the Lagrangian is in general given by 

1 
2 


PIy — hlaLi^NRoL + {Ms)paSLpNRa + ^ a^NRa^Ry + {Ps)pxSLpSix] + H.C. , 

( 6 ) 


where Ms is a Dirac mass term and fiR^s are Majorana mass terms. After EWSB, 
the Lagrangian (|^ gives rise to the following neutrino mass matrix in the flavor basis 
NRa, {Slp)^}- 

/ 0 Mn 0 \ 


M, = 


V 


0 

Ml 

0 


Md 

hi? 

Ms 


0 

MJ 




0 

M'^ 


D 


■Md 

■Mn 


(7) 


which has a form similar to the Type-I seesaw matrix (|^, with Md = (Mr),0) and 

, . / Mj ' 

Mn — 


1. Here we have not considered the dimension-4 lepton-number 
Ms fis ' 

breaking term which appears, for instance, in linear seesaw models 


134 


137 


since the mass matrix in presence of this term can always be rotated to the form given 
by ([^ [138 . Also observe that the inverse seesaw model discussed originally in 31 


set 


the RH neutrino Majorana mass /ij? = 0 in ([^. At the tree-level, the light neutrino 
mass is directly proportional to the Majorana mass term fis for \\fis\\ ^ ll^sll- 


M, = MrM-^iisM-^ Ml + 0{^l%), 


( 8 ) 


whereas at the one-loop level, there is an additional contribution proportional to 
yiR 41,139 , arising from standard electroweak radiative corrections 36 . The smallness 


of h'R,s is ‘technically natural’ in the’t Hooft sense 140 , i.e. in the limit of yiR^s 0, 


lepton number symmetry is restored and the light neutrinos vri are massless to all orders 
in perturbation theory, as in the SM. 
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Figure 1. Contours of the ratio of the average decay width of the heavy neutrinos to 
their mass splitting in the inverse seesaw model. The LNV signal will be unobservable 
in the shaded region with > 10. 


The freedom provided by the small LNV parameter /i^ in ([^ is the key feature of 
the inverse seesaw mechanism, allowing us to £t the light neutrino data for any value 
of light-heavy neutrino mixing, without introducing any hne-tuning or cancellations in 
the light neutrino mass matrix ® 141,142 . In essence, the magnitude of the neutrino 


mass becomes decoupled from the heavy neutrino mass, thus allowing for a large mixing 






10 


-2 


' 1 keV 
hs 


(9) 


The heavy neutrinos Nr and Sr have opposite CP parities and form a quasi-Dirac state 
with relative mass splitting of the order k = fis/Ms- All LNV processes are usually 
suppressed by this small mass splitting. For instance, in the one-generation case, the 
light neutrino mass in ([^ can be conveniently expressed as M,^ ~ KVi^Mu, in contrast 
with VinMo in the Type-I seesaw case [cf. Q]. It should be noted here that the 
approximately L-conserving models with quasi-degenerate heavy Majorana neutrinos 
could provide a natural framework 143 146 for realizing the mechanism of resonant 
leptogenesis 147-149 , where the leptonic CP asymmetry is resonantly enhanced when 
the mass splitting AMjy is of the same order as the decay width T^r. 

As for the LNV signature at colliders, in a natural seesaw scenario with approximate 
lepton number conservation, the LNV amplitude for the on-shell production of heavy 
neutrinos at average four-momentum squared s = -|- M%^)/2 can be written as 

2 AM R[ 


A 


LNV 




£N 


AM^+n 


+ o 


N 


/AMn 


M, 


N 


( 10 ) 


for AMjsi < r^, i.e. for small mass difference AMjq = \ between the heavy 

neutrinos compared to their average decay width T^r = (TjVi + T n 2 )I‘^- Thus, the LNV 
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Figure 2. Absolute eigenvalues of the neutrino mass matrix 0 as a function of the 
LNV mass term /is in the one-generation case with yn = 0 (left panel) and yn = ys 
(right panel). Here we have chosen = 1 GeV and Ms = 100 GeV for illustration. 
The horizontal dotted line shows the light neutrino mass scale required to satisfy the 
atmospheric neutrino data. 


amplitude in ( 10 ) will be suppressed by the small mass splitting, except for the case 
AMtv — Tat when it can be resonantly enhanced 147,150 . This suppression for the 


inverse seesaw case is illustrated in Figure where we show the contours of Fjv/AMjv 
for different values of the inverse seesaw mass parameters Mn and yis- As we can see 
from the plot, to observe LNV with M^r > 0(100) GeV, one needs a large /i 5 , and 
therefore, small [cf. ([^], which will suppress the LNV signal. Its implications for 

collider searches of these scenarios will be discussed in Section 12.31 

By virtue of the last equivalence in ([^, the inverse seesaw mechanism can be 
regarded as a variation of the Type-I seesaw mechanism. Therefore, it is instructive 
to study the transition between the neutrino mass formulas given by 0 and For 
illustration purposes, let us take a simplified version of ([^ for a single generation case. 
The neutrino mass spectrum for this scenario is shown in Figure as a function of the 
LNV parameter fis- As an example, we have chosen the Dirac mass term Mu = 1 
GeV and the heavy neutrino mass term Ms = 100 GeV such that Vim — Ms,/Ms = 
10 “^ is consistent with the current upper limit set by the electroweak precision data 
(EWPD) [^. First we consider the original inverse seesaw model with /i/j = 0 in ([^. 
As shown in Figure]^ (left panel), successful light neutrino mass generation occurs in this 
case only for /i ~ 1 keV [cf. ([^]. As evident from (|^, for jns <C M^, the lightest neutrino 
mass is proportional to whereas the two heavy neutrinos form a quasi-degenerate 
Dirac pair with masses Ms ± iis- For h-s ^ Ms, the heaviest neutrino mass becomes 
equal to p .5 and the lighter ones form a quasi-degenerate Dirac pair with mass of order 
Ms). For 7 ^ 0, the situation remains unchanged for the general case ^r^s ^ Ms, as 
shown in Figure]^ (right panel). However, for hr^s 3> Ms, we recover the Type-I seesaw 
with the lightest neutrino mass given by —M/^/jjiR [cf. (|^], whereas the heavier ones 
form a Majorana pair with masses equal to 
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Various laboratory searches have put stringent constraints on sterile neutrino mixing 
with active ones in a wide mass range of Mn from eV - TeV. For Mn values well 
below 1 MeV, as e.g. in eV-seesaw models EUl. the sterile neutrinos can be probed 
by neutrino-oscillation experiments. Assuming all sterile neutrinos to be of the same 
order of magnitude, current data rule out 1 neV < Mjsf < 1 eV 151 -154 . For 


10 eV < M]sf < 1 MeV, the mixing of sterile neutrinos with electron neutrino has 
been constrained by searches for Qufdfd and precision measurements of /d-decay energy 
spectra. For 1 MeV < Mjv ^ 1 GeV, the mixing with both electron and muon neutrinos 
have been constrained by peak searches in leptonic decays of pions and kaons. Sterile 
neutrino mixing with all neutrino flavors in the MeV-GeV mass range has also been 
searched for through their decay products in beam dump experiments. Upper limits 
on the active-sterile neutrino mixing elements have also been derived from cosmological 


bounds on sterile neutrino lifetimes as required for the success of BBN 155 157 


Here we summarize the current state-of-the-art sterile neutrino searches in the mass 
range 100 MeV < < 500 GeV, as relevant for collider experiments at the energy 

frontier and other planned experiments at the intensity frontier. Figures |3]j^ show the 
current constraints and some future projections on sterile neutrino mixing with the 
electron, muon and tau neutrinos, respectively. In these plots, the (gray) contour labeled 
‘BBN’ corresponds to a heavy neutrino lifetime > 1 sec, which is disfavored by BBN 
constraints 155-157 . The (brown) line labeled ‘Seesaw’ shows the scale of mixing as 
expected in the canonical seesaw [cf. d^]. We should remember that both these limits 


may get substantially modified in presence of more than two heavy neutrinos 
Other limits shown in Figures [3]j^ are explained below. 


133 


2.2.1. Neutrinoless Double Beta Decay The contributions of heavy Majorana neutrinos 
Njia to 0iyj3l3 amplitude is described by the standard neutrino exchange diagram between 
two /9-decaying neutrons, via a non-zero admixture of a weak eigenstate parametrized 
by the mixing element VeN^- The half-life is given by 


T'Ou 

1/2 


= A 


"V TJ. 


ei^i 


+ m. 


N 




rr {p"> + Ml 


( 11 ) 


where A = and (p^) = mpme\JCC^/. Here nrip and nrie are the 

proton and electron masses respectively, is the phase-space factor, is the nucleon 
axial-vector coupling constant, Uei is the light neutrino mixing matrix, ruj and Ma 
are respectively the light and heavy neutrino mass eigenvalues, and are 

the corresponding nuclear matrix elements (NMEs). These NMEs are conventionally 
calculated for the limiting cases pp (light) and S> pr (heavy), pp being 

the characteristic momentum transferred via the virtual neutrino, which is ~ 200 MeV 
corresponding to the mean nucleon momentum of Fermi motion in a nucleus. However, 
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jq-12- ^^^ ~ ■ .~~. ' . 

0.1 1 10 100 

Mn (GeV) 


Figure 3. Limits on the mixing between the electron neutrino and a single heavy 
neutrino in the mass range 100 MeV - 500 GeV. For details, see text. 


the interpolating formula ( pT| allows us to calculate the 0z//3/3 half-life for arbitrary 
heavy neutrino masses using the NMEs (light) and (heavy) (T^[T^ . 


Using the combined 90% C.L. limit on Oiy/d/d half-life > 3 x 10^® yr from 

GERDA-I-Heidelberg-Moscow experiment 160 , we derive from the second term in ( [IT| ) 
upper limits on |UeAfP as a function of a generic heavy neutrino mass Our results 
are shown in Figure where the shaded (orange) region between the solid and dashed 
lines, labeled ‘GERDA’, shows the uncertainty due to NMEs 159,161 . Here we have 
used the recently re-evaluated phase-space factors 162 and the NMEs from a recent 


calculation within the quasi-particle random phase approximation (QRPA) [159 163 


Similar limits are obtained using the half-life limit ^ 2.6 x 10^® yr from 


KamLAND-Zen experiment 164 165 and the corresponding QRPA NMEs 159 


From Figure]^ it seems that the Oi^/d/d constraints are very severe, thus shadowing 
the future prospects of observing LNV in other processes involving the electron channel. 
However, one must keep in mind that the OufdlS limits may be signihcantly weakened 
in certain cases when a cancellation between different terms in ( pH| may happen 166 , 
e.g. due to the presence of Majorana CP phases. In general, the Majorana nature of 


neutrinos does not guarantee an observable rate in all models 167 . Also, in the 


inverse seesaw scenario with pseudo-Dirac heavy neutrinos, the Ov/d/d limits are usually 
diluted by the small LNV term k = fis/Ms- Therefore, it is still important to include 
the electron channel while performing an independent direct search for heavy neutrinos 
at colliders. 
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Figure 4. Limits on the mixing between the muon neutrino and a single heavy 
neutrino in the mass range 100 MeV - 500 GeV. For details, see text. 


2.2.2. Peak Searches in Meson Decays Peak searches in weak decays of heavy leptons 
and mesons are powerful probes of heavy neutrino mixing with all lepton flavors. The 
most promising are the two-body decays of electrically charged mesons into leptons and 
neutrinos: ^ i^N 168 170, whose branching ratio is proportional to the mixing 


iV^ATp. Thus, for a non-zero mixing and in the meson’s rest frame, one expects the 
lepton spectrum to show a second monochromatic line at 


Ef = 


Mj. -b m'j 


Ml 


2M 


X 


( 12 ) 


apart from the usual peak due to the active neutrino un. For sterile neutrinos heavier 
than the charged lepton, the helicity suppression factor inherent in leptonic decay rate is 


weakened by a factor Ml/mj 169 due to which the sensitivity on | increases with 
Mat till the phase space becomes relevant. Peak searches have been performed in the 
channels TT —>• eN 171-175 , vr —)■ /xiV [176 -180 , K —>■ eN 181 and iP —)■ 181-185 


The current 90% C.L. limits on \Vm\^ (for t = e,/i) derived from these searches are 
shown in Figures]^ andlabeled as ‘X —)■ (with X = ti,K and i = e, fx). The limit 
from 71 —)■ fiN is not shown here, since it is only applicable in the mass range 1 MeV 
^ M]^ ^ 30 MeV. 

The peak searches could in principle be extended to higher masses with heavier 


meson/baryon decays 186 -188. For instance, the Belle experiment 189 used the decay 


mode B —)■ XiN followed by iV —)■ f'vr (with £ = e, /i) in a data sample of 772 million 
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Figure 5. Limits on the mixing between the tan neutrino and a single heavy neutrino 
in the mass range 100 MeV - 500 GeV. For details, see text. 


BB pairs coming from T(4s) resonance to place 90% C.L. limits on |14AfP and in 

the heavy neutrino mass range 500 MeV to 5 GeV, as shown in Figures and labeled 
as ‘Belle’. 

Limits on the mixing parameter can also be set from the 3-body decay of muons, 


where a sterile neutrino contribution may distort the spectrum of Michel electrons 169 
In case of r-leptons, the 2-body decays into hadrons r —?• NX are promising. If the 
hadronic system X hadronizes into charged pions or kaons, then its mass and energy 
can be reconstructed at high precision. Using future i?-factories with a large dataset of 
r decays like t~ —)■ , stringent limits on the mixing parameter lU-AfP can be 

placed 190 , as shown in Figure]^ where the (red, solid and dashed) contours labeled 
B-factory are the conservative and optimistic projected limits at 90% C.L. from ~ 10 
million r-decays. 

We should note here that the bounds from peak searches are very robust since they 
use only the kinematic features and minimal assumptions regarding the decay modes 
of the heavy neutrino. Moreover, since the heavy neutrino is assumed to be produced 
on-shell, these limits are valid irrespective of whether the heavy neutrino is a Majorana 
or Dirac particle. 


2.2.3. Beam Dump Experiments Another way to constrain the sterile neutrinos is via 
searches of their decay products. The sterile neutrinos are unstable due to their mixing 
with active neutrinos, and their decay rate is proportional to the mixing parameter 
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iV^ATp. Thus, if kinematically allowed, they can be produced in semileptonic meson 
decays, and then subsequently decay into visible final states such as charged leptons, 
pions and kaons. These visible products can be searched for in beam dump experiments 
by placing the detector some distance away from the production site. The current 90% 


C.L. limits from some of these beam dump experiments, such as PS191 191 , NAS 192 
CHARM flgsIpS) , IHEP-JINR [^, BEBC 


197 , FMMF 198 , NuTeV 199 and 


NOMAD 200 
CHARM 


are 


shown in Figures [3]|^ It may be noted here that the PS191 [191 and 


193 limits shown here assume that the sterile neutrinos interact only via CC. 


Including the NC interactions requires a reanalysis of the data which, in the context of 
i/SM, gives twice stronger limits 1201]. 


The proposed LBNE (now DUNE) experiment with a near detector could probe 
much smaller values of the mixing |202| , as illustrated in Figures |3]|^ by the projected 
limits labeled ‘LBNE’. Here the heavy neutrinos are produced in charmed meson decays, 
and a near detector length of 30 m with a ~ 5 year exposure of 5 x 10^^ protons on 
target is assumed. We have obtained the limits on individual mixing parameters from 
the corresponding limit on the sum of mixing iWvP assuming a normal hierarchy of 
light neutrinos. 


2.2.4- R<ife LNV Decays of Mesons For heavy Majorana neutrinos, one could also 
look for rare LNV decays of mesons: —)■ N —)■ which are forbidden 

in the SM. Searches for such decay modes have been performed in many experiments, 
such as CLEO, Belle, BaBar and LHCb [^. The most stringent constraints come from 

shown in Figures and Here a realistic detector size of 


K~^ —)■ mode 203 


as 


10 m has been assumed. The corresponding limits from D and B meson decays 204 206 


are found to be weaker than the existing limits in the relevant mass region, and therefore, 
are not shown here, except the recent LHCb bound 207 in Figure]^ which was obtained 
using the B~ —)■ decay mode with 3 fb“^ of integrated luminosity collected at 

/s = 7 and 8 TeV LHC. 

The bounds discussed in Sections [2.2.3l and l2.2.4l are less robust than those discussed 


in Section |2.2.2| because they are weakened, or even completely evaded, if the sterile 
neutrinos have other dominant decay modes into invisible particles. If the sterile 
neutrino decay length is shorter than the detector size, the number of signal events is 
suppressed by lUwP- On the other hand, if their decay length is larger than the detector 
size, the sterile neutrinos decay mostly outside the detector and the number of events is 
further suppressed by | This limitation could be overcome by increasing the flux of 


initial hadrons, e.g. in proposed fixed-target experiments such as SHiP 208 using high- 


intensity proton beams at the CERN SPS. The huge background due to multiparticle 
production inherent in hadron scatterings can be absorbed by adopting appropriate 
beam-dump techniques, thus allowing the sterile neutrinos to freely propagate into a 
decay volume. This experiment will improve the mixing sensitivity by up to four orders 
of magnitude 209 , as shown by the projected limits in Figures [3]j^ 
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2.2.5. Z-decays For < Mz, using the possible production of heavy neutrinos in 
the Z-boson decay Z —>■ or Z —)■ 210| , and its subsequent CC and NC 

decays, 95% C.L. limits on the mixing parameters |VfArp were obtained by L3 211 
and DELPHI 212 collaborations from a reanalysis of the LEP data. These limits 


are 


shown by the contours labeled ‘L3’ (pink, dashed) in Figures and and by the 
contours labeled ‘DELPHI’ (dark green, dashed) in Figures |3]j^ 

A future high-luminosity Z-factory, such as the proposed FCC-ee experiment, 
will dramatically improve the sensitivity down to ~ 10“^^ for mixing with all 

neutrino flavors and covering a large phase space for heavy neutrino masses in the 10-80 
GeV range 213,214. This is shown in Figures [sll^ where the limits on individual 


mixing parameters are derived from the corresponding limit on the sum of mixing 
IV^TVP 213 assuming a normal hierarchy of light neutrinos. We have also assumed 
10^^ Z-boson decays occurring between 10-100 cm from the interaction point. Increasing 
the number of Z-bosons and/or the range of decay length could further enhance these 
sensitivity limits, eventually reaching the theoretical expectation from the canonical 
seesaw formula, as shown by the (brown) dashed line labeled ‘Seesaw’ in Figures |3]j^ 


2.2.6. Electroweak Precision Tests Due to their mixing with active neutrinos, heavy 
neutrinos can affect various EW precision observables, such as the Z invisible decay 
width and electroweak parameters in the SM 215 -217 . The same mixing effects also 


show up in non-unitarity of the leptonic mixing matrix 66 - 68 and the violation of lepton 
universality in leptonic and semileptonic decays of pseudoscalar mesons |218 220 . Using 


global hts to the EWPD, stringent model-independent constraints on have been 

derived 69 73 . The current 90% C.L. limits are are shown in Figures |3]j^ (brown, dotted 
contours labeled ‘EWPD’). These limits are independent of the heavy neutrino mass for 
M^v > Mz, and there is a mild mass dependence for lower Mjv values. Here we have 
only included the electroweak precision observables and lepton universality observables 
in the ht. The LFV observables are not included here, since they are more sensitive to 
the details of the Yukawa structure in the underlying model. 


2.3. Direct Collider Searches 

Heavy neutrinos with masses of the order of electroweak scale can be directly produced 
on-shell at colliders. Such a direct search was performed in annihilation at 


LEP 221,222 , assuming a single heavy neutrino production via its mixing with active 
neutrinos: —)■ Nfu, followed by its decay via NC or CC interaction to the SM 

W, Z or Higgs {H) boson: N — )■ iW, ^mZ, Concentrating on the decay channel 

N —)■ eW with W —)■ jets, which would lead to a single isolated electron plus hadronic 
jets, the L3 collaboration put a 95% C.L. upper limit on the mixing parameter |I4v| 


2 ; 


m 


a heavy neutrino mass range between 80 and 205 GeV [222 , as shown by the (red, solid) 
contour labeled ‘L3 IP in Figure This search was mainly limited by the maximum 
center-of-mass energy y/s = 208 GeV at LEP. Future lepton colliders can signihcantly 
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improve the sensitivity in this mass region, as illustrated in Figure by the projected 
limit labeled ‘ILC’, which is obtained assuming a \/s = 500 GeV ILC with luminosity 
of 500 fb"^ 


223 


In the context of hadron colliders, a Majorana heavy neutrino leads to the smoking 
gun lepton-number violating signature of same-sign dilepton plus jets with no missing 
transverse energy: pp {pp) —)■ W* —)■ Ni^ —>• 36,1^203,224-229 . An inclusive 


search for new physics with same-sign dilepton signals was hrst performed in pp collisions 
at the Tevatron 230 . After the inauguration of the LHC era, the CMS and ATLAS 


collaborations have performed direct exclusive searches for the on-shell production of 
heavy neutrinos above the Z-threshold. The previous searches with 4.7 fb~^ data at 
y/s = 7 TeV LHC set 95% C.L. limits on ^ 10“^ — 10“^ (with i = e, p) for heavy 

neutrino masses up to 300 GeV 231, 232| . More recently, these limits were extended 
for masses up to 500 GeV with 20 fb“^ data at y/s = 8 TeV 233,234 , and are shown 
in Figures]^ (ATLAS) and(ATLAS and CMS). For Mat ~ 100 GeV, the direct limits 
in the muon sector are comparable to the indirect limits on |V)(Arp ^ 10“^ imposed by 
the electroweak precision data 69,72 and LHC Higgs data [235 , 236| . With the run-H 
phase of the LHC starting later this year with more energy and higher luminosity, the 
direct search limits could be extended for heavy neutrino masses up to a TeV or so. 
Also note that the LFV processes put stringent constraints on the product 
(with i 7 ^ i') [237 -240 , but do not restrict the individual mixing parameters |V^Arp in a 
model-independent way. So the direct searches provide a complementary way to probe 
the light-heavy neutrino mixing in the seesaw paradigm. 

All the direct searches at the LHC so far have only considered the simplest 


production process for heavy neutrinos through an s-channel hF-exchange 36,150,203 
225-229 , as shown in Figure 


However, there exists another collinearly enhanced 
electroweak production mode involving t-channel exchange of photons: pp —)■ W*^* —>■ 
Ni^jj 241 , cf. Figure 0 which gives a dominant contribution to the heavy neutrino 
production cross section for higher Mn values. This is mainly because of the fact that 
with increasing heavy neutrino mass, the production cross section for the t-channel 
process drops at a rate slower than that of the s-channel process, though the exact 
cross-over point depends crucially on the selection cut for the pt of the additional jet 
associated with the virtuality of the t-channel photon 241,242 . The photon-mediated 
process pp —)■ W*'^* —)■ Ni^jj has two contributions: an inelastic part with t-channel 
virtual photon and an elastic part py —)■ Ni^j with a real photon emission from one of 
the protons. The elastic part is calculated using an effective photon structure function 
for the proton 243,244 , whereas the inelastic part is computed for a non-zero minimum 
Pt of the jet associated with the virtuality of the photon [241 . A comparison of the 
production cross sections for the s-channel Drell-Yan (DY) process pp —)■ W* —)■ N 
and the photon initiated processes is made in Figure for a representative value of 
p^min ~ 20 GeV which, alongwith a jet separation cut NRP > 0.4, is sufficient to ensure 
that the photon-mediated processes are collinear safe. We hnd that the total photon- 
initiated contribution becomes dominant over the DY cross section for Mat > 600 GeV. 
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Figure 6. Feynman diagram for heavy neutrino production at the LHC via the s- 
channel Drell-Yan process. 



Figure 7. Feynman diagrams for heavy neutrino production at the LHC via the 


t-channel photon-mediated processes 241 


Here, we have not included the QCD corrections, which could further lower the cross¬ 
over point to the level presented in 241| , but this requires a more careful analysis and 
will be presented elsewhere. Note that the numerical results shown in Figure are 
slightly different from those presented in 242 , which can be mainly attributed to the 
choice of regulator used to treat the collinear behavior. 

In any case, including the collinear enhancement effect could further enhance the 
heavy neutrino signal sensitivity at the next run of the LHC 241 . As an illustration. 


we have shown in Figures and projected conservative limits with 300 fb ^ data at 
^/s = 14 TeV (blue, dashed contours labeled ‘LHC 14’), assuming that the cross-section 
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Figure 8. Comparison of the cross sections for heavy neutrino production at ^/s = 14 
TeV LHC via the s-channel (Figure]^ and t-channel (Figurediagrams. 


limits are at least as good as the existing ones at a/s = 8 TeV, as reported in 234 


The direct collider limits for < 100 GeV are not likely to improve signihcantly with 
higher collision energy, due to the increased pile-up effects, thus obfuscating the low-pr 
leptons produced by the decay of a low-mass heavy neutrino. Instead, a displaced vertex 
search might be useful to probe the low-mass range between 3-80 GeV for mixing values 
10 -' 10 -^ 1245 . 


For heavy Dirac neutrinos as predicted in theories with approximate L-conservation 
[cf. (|^], the same-sign dilepton signal is suppressed. In this case, the golden channel 
is the trilepton channel; pp —)■ W* —)■ —)■ 246- 


252 


Using this 


trilepton mode and also taking into account the infrared enhancement effects 241 


direct limits on the mixing of heavy Dirac neutrinos with electron and muon neutrinos 


were obtained 251 by analyzing the tri-lepton data from ^/s = 8 TeV LHG 253 


Finally, we note that there exist no direct collider searches for heavy neutrinos 
involving tan-lepton hnal states. This is mainly due to the experimental challenges of r 
reconstruction at a hadron collider. The situation is expected to improve in future with 
better r-tagging algorithms and/or in cleaner environments of a lepton collider. 


3. Heavy Triplets at Colliders 

Unlike the minimal Type-I seesaw messengers which, being SM gauge singlets, can only 
communicate with the SM sector through their mixing with the active neutrinos, the 
Type-II and III seesaw messengers are SU{2)l triplet scalar (A’*"’'', A"*", A°) and fermion 
helds respectively, and hence, can be directly produced at the LHG via 


their gauge interactions. For Type-II seesaw 15 -19 , the smoking gun signal would be 


the detection of a doubly-charged scalar with LNV interactions. For this scenario, the 
most relevant production channels at the LHG are pp —)■ Z* A+’*’A , A+A“, 

pp and pp ^ W* ^ A=^=^A=f, A=^=^IU^ |^[254[j^ . The 
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doubly-charged scalar boson has the following possible decay channels: 

and A^A^, if kinematically allowed. For the triplet VEV ua 0.1 MeV, the 
doubly charged Higgs couplings to is suppressed and for a nearly degenerate triplet 
mass spectrum, the dominant decay mode of A^^ is same-sign dileptons 260,265 . In 


this case, the current 95% C.L. experimental lower bound on the doubly-charged triplet 


scalar mass is between 490-550 GeV, depending on the final lepton flavor 266-268 


For ua ^ 0.1 MeV, the Yukawa couplings of the A fields to leptons is suppressed and 
A^^ dominantly decays to same-sign dibosons, in which case the collider limits are 
significantly weaker 269 -273 . These mass bounds could be significantly improved in 


the upcoming run II phase of the LHC. A future lepton collider such as the ILC will 
offer an almost background-free environment for the doubly-charged scalar detection, 
if kinematically allowed, through the single production process e’''e“ —)■ e^£^A^^ 274 


or pair-production process e+e —)■ A’'''^A 275 . We note here that a relatively 


light charged scalar will affect the loop-induced decays of the SM Higgs boson h —?• 
and h —)■ Zj 


260,276-281 


278-282 . In fact, for a given enhancement in these 


77 

decay rates over the SM prediction, one can set an upper limit on the charged scalar 
mass in the minimal Type-H seesaw model using vacuum stability and perturbative 
arguments 279 . 

the relevant production mechanisms at the LHC are 


20 


For Type-HI seesaw 

pp Z* j-i* S+S", pp^ Z* ^ and pp 

The neutral component has the same decay modes as the heavy sterile neutrino 


246 283-285 


discussed in Section 2.2, i.e. —)■ h'i{i'e)Z, , whereas the charged fermion 

has the LNV decay modes V(,{V()W^,C^Z,C^H. The current 95% C.L. experimental 


lower bound on the fermion triplet mass is 245 GeV 286|287 . This model has promising 


discovery prospects at the ILC as well 288,289 


4. Extended Gauge Sectors 

Within the SM gauge group, there is no explanation for the origin of the Majorana 
masses of the seesaw messengers. This problem can be solved by extending the SM gauge 
group so that the Majorana mass can be associated with the spontaneous breaking of the 
extra gauge symmetry. The simplest of such additional gauge symmetries is U{1)b-l 
whose breaking could set the mass scale of the RH neutrinos in the Type-I seesaw. A 
natural embedding of the RH Majorana neutrinos as well as the scalar triplets can be 
found in the LRSM based on the gauge group SU{2)i x SU{2)b x U{1)b-l- Therefore, 
we will focus on these two additional gauge symmetries in the following. 


%i. Additional U{1) 

The simplest extension of the SM gauge group to explain the heavy Majorana neutrino 
mass in (j^ is the inclusion of an additional U{1) gauge group along with an associated 
Z' gauge boson. This symmetry can be added by hand to the SM but it could also 
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Mn (TeV) 


Figure 9. Left: Feynman diagram for heavy neutrino production via the Z' resonance. 
Right: Heavy neutrino decay length as a function of its mass Mn and mixing Vin 
(solid blue contours). The dashed red contours denote Br(/i —)• ey), with the shaded 
red region on top excluded by the current MEG limit [294| . The gray shaded band 
highlights the parameter range where light neutrino mass scales between and 

0.3 eV are generated within the canonical type-I seesaw mechanism. 


naturally arise from a UV-complete theory such as the Pati-Salam model, S'O(IO) or 
Eq GUTs. For a review of various Z' models, see e.g. 290 . The mass and couplings of 


the Z’ boson are strongly constrained by EWPD. Constraints from lepton universality 
at the Z peak puts a lower limit of Mz' > C){1) TeV 291 , whereas direct searches at 
the LHC exclude Mz' below about 2 TeV [^. Similarly, the mixing angle between Z' 
and the SM Z is limited to be less than 0(10“^). 

With regard to heavy neutrinos, the main phenomenological advantage is the 
possibility that the heavy neutrinos are charged under the additional U{1). This would 
provide a new and potentially strong production channel at colliders. A Feynman 
diagram for resonant heavy neutrino production at the LHC via such a Z' portal with 

shown in Figure (left panel). An 


a hnal state of two leptons and four jets 292 


IS 


important point to note is that the total cross section of this process is independent of 
the mixing strength Vin- This process could therefore be observed at the LHC as long 
as the total decay width of the heavy neutrino is large enough such that it decays within 
the detector. The relevant parameter range for this to occur is shown in Figure]^ (right 
panel) [293 . Even the canonical Type-I seesaw with small mixing |VfAr| ^ 10“® for TeV- 
scale Mtv can be potentially probed in this case, possibly through displaced vertices. 
This also includes the potential to observe LFV signatures, despite the unobservably 
small LFV rates for low-energy processes such as /i —)■ ey, as they are strongly suppressed 
by such a small mixing. 
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4 . 2 . Left-Right Symmetric Models 

The minimal LRSM which extends the SM gauge symmetry to SU{2)l x S'f/(2)/j x 
U( 1 )b-l [^[^[ 2 ^ provides a simple UV-complete seesaw model, where the key 
ingredients of seesaw, i.e. the RH neutrinos and their Majorana masses, appear 
naturally. The presence of RH neutrinos is a necessary ingredient for the restoration of 
Left-Right symmetry and is also required by anomaly cancellation, whereas the seesaw 
scale is identihed as the breaking of the S'f/(2)^ symmetry. It is worth noting that in 
the presence of three RH neutrinos, the (R —L)-symmetry which was a global symmetry 
in the SM becomes a gauge symmetry in the LRSM, as the gauge anomalies cancel by 
satisfying the condition Tr(R — L)^ = 0. Moreover, the electric charge formula takes a 
form similar to the Gell-Mann - Nishijima relation: Q = hr + hn + [B - L)/2, where 
Jsi and are the third components of isospin under SU{2)l and SU{2)r respectively, 
and the SM hypercharge can now be understood as 1^/2 = hR + {B-L)/2 |^[^ . 

In the LRSM, leptons are assigned to the multiplets = (z/^, t )r and = (A^^, €)r 
(where i = e, fi,T is the generation index) with the quantum numbers ( 2 , 1 , — 1) and 
( 1 , 2 , —1) respectively under SU{ 2 ) x SU{ 2 ) x U{1)b-l- The Higgs sector of the minimal 
LRSM contains a bidoublet 0 with quantum numbers (2,2,0) and two triplets Al,r 
with quantum numbers ( 3 , 1 , 2 ) and ( 1 , 3 , 2 ) respectively. The VEV vr of the neutral 
component of Ar breaks the gauge symmetry SU{2)r x U{1)b-l to U{1)y and gives 
masses to the RH gauge bosons Wr, Zr boson and the RH neutrinos Nr. The VEVs 
(k, k') of the neutral components of the bidoublet 0 break the SM symmetry and are 
therefore of the order of the electroweak scale. 

The LRSM Lagrangian relevant for the neutrino mass is given by 


-Cy = hLcfR + hL(pR + fLL'^Cia 2 AiL + /rR^ Cia 2 ARR + H.c. , 


(13) 


where cj) = a 2 (j)*o '2 and h, h, fi^R are 3x3 complex Yukawa couplings. After symmetry 


breaking, (13) leads to the Dirac mass matrix Mr = hn + hn' and the Majorana mass 


matrices Mr = fiVn and Mr = JrVr for the light and heavy neutrinos respectively, 
where ul is the VEV of the neutral component of This leads to the neutrino mass 
matrix 


M,. = 


Ml 
Ml 



(14) 


as compared to the Type-I case given by ([^. In the usual seesaw approximation 


IlM^ill <C ||Vf^||, diagonalizing (14) leads to the light neutrino mass matrix of the form 

M. = Mr, - MnM^^Ml , (15) 


D ) 


where the second term on the RHS is the Type-I seesaw contribution which is inversely 
proportional to vr, whereas the first one is the Type-H seesaw contribution which is 
directly proportional to vl- It should be noted here that in the minimal LRSM, if charge 
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conjugation is the discrete L-R symmetry, Mr, = Mj, and Ml = {vlIvr)Mr. In this 
case, the Dirac Yukawa couplings are generically constrained by the light and heavy 
neutrino mass and mixing parameters 298 . However, there are exceptions, e.g. in H 4 
symmetry-based models, where 


( 


Ml oc 


V 


a; 


UJ 

2 


1 


bJ 

1 




/ 


(16) 


with = 1 and Mat = mArla. For such symmetry-based models. My vanishes identically 
for ua = 0, independently of the size of the Dirac Yukawa couplings]^ Likewise, in 


versions of LRSM where parity and SU{2)ji gauge symmetry scales are decoupled 299 


the Dirac Yukawa couplings could be sizable 132 , while being consistent with the light 
neutrino data. 

Since processes induced by the RH currents and particles have not been observed so 
far, vr has to be sufficiently large. In particular, hadronic flavor changing neutral current 
effects restrict Mwj^ ^ 3 TeV |300 303 , assuming that the SU{2)r gauge coupling gR 
has the same strength as the SU{2)l gauge coupling ql- Direct search limits from the 
a/s = 7 and 8 TeV LHC data put similar constraints on M^^, depending on the heavy 

This translates into a lower limit on vr = ^/2Mwp./Q r ^ 6.5 


neutrino mass 


TeV. On the other 


304 


305 


land, for the the left-triplet VEV vr, the electroweak p-parameter 


constraints set an upper limit on < 5 GeV 306 


Due to the presence of RH gauge interactions, the LRSM gives rise to a number of 
new contributions to both LNV and LFV processes; see e.g. 307 309| . In particular, 
there are several diagrams that contribute to the amplitude: (i) standard light 

neutrino exchange with mass helicity flip 310 311| , (ii) RH neutrino and RH gauge 
boson exchange 12 ,[^, (iii) RH Higgs triplet exchange [312 


and 


IV 


mixed LH-RH 


contributions 78 309,313- 316j . The latter depend on the size of the left-right neutrino 
mixing Vm — MrM]^^. In the Type-H seesaw dominance 317,318 , neglecting the 


mixed left-right as well as the canonical light-neutrino contributions, the Offd/d half-life 
due to purely RH currents can be written as 


'~rOu 

1/2 


= A 


TUn 


Ml 


w 


M, 


Wr 


T/2 

M, 


( 17 ) 


where 14* is the mixing matrix for the RH neutrinos Ni with mass eigenvalues Mj and 
A is defined below ( 11 ). Using 0 and the current experimental limits on T^, lower 


limits on the RH gauge boson and RH neutrino masses can be derived 161,307 . This 


is illustrated in Figure 10, where the excluded areas from Ouldfd searches are shown. 
Similarly, low-energy LFV processes such /i —>■ ey and n ^ die can be drastically 
enhanced in the LRSM, with a host of new contributions 319 . This is also illustrated 


in Figure 10 using maximal e/i flavor mixing of the heavy neutrinos 307 


^ If parity is the discrete L-R symmetry which implies Mr = M/ , a similar construction can be made 
by interchanging the second and third columns of Mr in (161. 
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Mw, (TeV) 


Figure 10. Left: Feynman diagrams contributing to the ‘smoking gun’ collider signal 
of LNV {(-^tpjj) in the LRSM through the production via SM kF(L) and heavy Wr, 
giving rise to 4 different contributions: RR, RL, LL, LR. Right: Comparison of LNV 
event rates via the RR diagram at the LHC and in Oi//3/3 experiments 307 . The solid 


blue contours give the signal signihcance of 5cr and 90% at the LHC with 14 TeV and 
C = 300 fb“^. The area denoted ‘LHC excL’ is excluded by current LHC searches 
in the electron channel 305 . The green contours show the sensitivity of current and 


future 0vj5f5 experiments, assuming dominant doubly-charged Higgs or heavy neutrino 
exchange and the red contours show the sensitivity of LFV processes as denoted. 


As for the LHC phenomenology, the presence of RH gauge interactions could lead 
to signihcant enhancement of the LFV/LNV signal. There are several contributions 
to the smoking gun LNV signal of same-dilepton plus two jets, as summarized in 
Figure 10 (left). Even if the left-right neutrino mixing is small, heavy RH neutrinos 
could be directly produced via s-channel Wr exchange and subsequently decay via the 
same Wr 224 . The potential to discover LFV and LNV at the LHC in this scenario 


has been analyzed in 307,320 - 323| . Figure [l0| compares the sensitivity of LNV searches 
at the LHC with the sensitivity of OujSjS experiments. The y/s = 14 TeV LHC might be 
ultimately able to probe RH gauge boson masses up to Mwj^ = 6 TeV 324 , whereas 


a futuristic ^/s = 80 (100) TeV hadron collider could probe up to = 26.6 (35.5) 


TeV 325 


In variations of low-scale LRSM with large left-right neutrino mixing 132 


there will be new contribution to the like-sign dilepton signal due to mixed RH-LH 
currents 326,327 , in addition to the purely RH and LH contributions. Note that the 


amplitude for the RR diagram in Figure(left) is independent of Vjv, and hence, does 
not probe the full seesaw matrix ( [I^ . On the other hand, the RL diagram is sensitive 
to the heavy-light mixing 326| , and in fact, the dominant channel over a fairly large 
range of model parameter space, as illustrated in Figure pTj Thus, a combination of the 
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Figure 11. Phase diagram illustrating the dominance of different channels, namely, 
LL, RL and RR as shown in Figure 10 (left panel) in different regions of LRSM 
parameter space. Herer we have chosen Mn = 100 GeV for illustration. The gray 
shaded region labeled ‘Seesaw’ is the scale of mixing parameter as expected in 

canonical seesaw [cf. &■ The brown shaded region labeled ‘EWPD’ is the typical 
range of mixing ruled out from electroweak precision data (see Section 2.2.6). 


diagrams shown in Fignrej^ (left) is essential to fnlly explore the seesaw mechanism at 
the LHC. In this context, it is also nsefnl to distinguish the RH gauge boson contributions 


to the collider signatures from the LH ones using different kinematic variables 3261328 


Recently, the CMS collaboration analyzed the a/s = 8 TeV LHC data with 19.7fb“^ 
resulting in the most stringent direct bounds on the RH gauge boson masses up to 

While the analysis hnds no signihcant departure from SM 


M, 


Wr 


= 3 TeV 305 


expectations, the LHC data exhibit an intriguing excess in ee production with a local 
signihcance of 2.8cr for a candidate Wr mass of 2.1 TeV; no excess is observed 

in the channel. The excess could be interpreted as a hint for Wr production, but 
with a smaller RH gauge coupling Qr/Ql — 0.6 329-332 . If the excess turns out to be 


statistically signihcant in future with more data and independent scrutiny from ATLAS, 


it might be an evidence for L-R symmetry with high-scale parity breaking 299 


5. Lepton Number Violation at the LHC and Leptogenesis 


The observed baryon asymmetry of the Universe 57 is far too large compared to the 


SM expectation. This is because in the SM, the necessary CP violation is too small and 
no hrst order phase transition can take place for the observed value of the Higgs mass. 
The mechanism of baryogenesis in the Universe thus requires physics beyond the SM. A 
large number of possible mechanisms to generate the observed baryon asymmetry have 
been proposed in the literature. A particularly interesting scenario that also connects to 


the neutrino mass mechanism is leptogenesis 333 . In its original formulation, the out- 
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of-equilibrium decay of the same heavy Majorana neutrinos responsible for the Type-I 
seesaw mechanism would create a lepton asymmetry, which is then reprocessed into 


a baryon asymmetry through [B + L)-violating EW sphaleron transitions 334 at or 
above the scale of EWSB. In ‘vanilla’ leptogenesis with hierarchical heavy neutrino 
335 , the neutrino oscillation data impose a /ower bound of M^r > 5 x 10^ GeV 

which is inaccessible to foreseeable 


masses 


336 338 


on the lightest RH neutrino mass 
laboratory experiments. Moreover, these high-scale thermal leptogenesis scenarios 
are in conflict with the upper limits on RH neutrino masses from EW naturalness 
arguments 83 86 , and in many supergravity models, from the gravitino overproduction 


bound 87 . These problems can be naturally avoided in the framework of resonant 


leptogenesis 147-149 , where the heavy-neutrino self-energy effects on the leptonic CP 


asymmetry become dominant 339,340 and get resonantly enhanced, even up to order 
1 147 148 , when two heavy Majorana neutrinos have a small mass difference comparable 


to their decay widths. This allows successful thermal leptogenesis with low seesaw scale 
accessible to laboratory experiments, while maintaining agreement with other theoretical 


and experimental constraints 121,128,341-343 


In addition, the observation of LNV at the LHC would have important consequences 
on the viability of general leptogenesis models. The issue of probing leptogenesis at the 


LHC has been studied in the context of the LRSM 344 346 but it is also possible 


to falsify a large class of high-scale leptogenesis scenarios if LNV was observed at the 


LHC 347,348 . The analysis of 347 focuses on the resonant LNV process pp —)■ l^l^jj 


involving generic intermediate particles as shown in Figure 12 (left). Such a class of 
diagrams is generated by general decompositions of the corresponding 9-dimensional 
short range Ouldfd operator [349 . As an example, the specific realization in the LRSM 
is shown in Figure [l0| (left). The minimal lepton asymmetry washout rate Tw induced 
by the process in Figure [I^ (left) is then related to the corresponding LHC cross section 
ctlhc as |M7 


logio ^ ^ 6.9 + 0.6 


Mx 

TeV 


1 + log 


10 


O'LHC 

liT’ 


( 18 ) 


where H is the Hubble parameter at the scale Mx, i.e. the mass of the resonance in 
Figure 12 (left). If Tw/H S> 1, the dilution of a primordial net lepton number density. 


understood to be produced by a leptogenesis mechanism at a higher scale, is highly 
effective and the lepton asymmetry would be washed out before it can be converted by 
sphaleron processes. This result is illustrated graphically in Figure 12 (right). Both 


(18) and Figure 12 (right) demonstrate that the observation of LNV at the LHC 
necessitates a very large lepton asymmetry washout. It would therefore rule out or 
strongly constrain leptogenesis scenarios above the scale Mx- Low scale scenarios, such 
as resonant leptogenesis discussed above where the lepton asymmetry is generated at 
scales lower than Mx, would not necessarily be constrained. 

The approach is not limited to resonant same sign dilepton production but can be 
applied to any other process with AL ^ 0, A{B — L) 7 ^ 0 and no missing energy at 
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Figure 12. Left: Generic diagrams for resonant same sign dilepton production 
pp —> CCqq at the LHC. The intermediate particles are generic scalar or vector 
bosons X, '> and a fermion Any two of the final state fermions p can be leptons, 
depending on the transformation properties of the intermediate particles under the SM 
gauge group. Right: Induced lepton number washout rate Vw/H at T = Mx as a 
function of the LNV scale Mx and the LHC cross section ctlhc (solid blue contours). 
The red dashed curves denote typically expected cross sections for gauge strength 
interactions. The shaded area is excluded by current LHC searches 


305 


the LHC, and potentially at other fntnre colliders. One example would be resonant 
pair-production of heavy particles, e.g. leptoquarks decaying to same-sign dilepton plus 


two jets 350 or singlet fermions decaying into a six fermion hnal state [cf. Figure]^. 
Searches for high-energy LNV processes can therefore have a signihcant impact on 
models of leptogenesis, and baryogenesis in general. 


6. Supersymmetric Models 

In spite of the null results for SUSY searches so far at the LHC, SUSY still remains 
one of the most popular BSM scenarios due to its many attractive features, such as 
the gauge coupling unihcation, radiative EWSB, amelioration of the gauge hierarchy 
problem, natural DM candidates, connection to gravity, and so on (for a review, see 
e.g. 351| ). However, the simplest supersymmetric extension of the SM, namely, the 
Minimal Supersymmetric SM (MSSM) with conserved i?-parity, does not accommodate 
non-zero neutrino masses. There are two common approaches for breaking the L- 
symmetry of MSSM to incorporate neutrino masses, and we briefly discuss their collider 
signatures in the following subsections. The possibility of connecting neutrino masses 

352H35^ ). 


to SUSY-breaking sector has also been investigated (see e.g. 
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Just as light neutrino masses can be generated in a seesaw extension of the SM, 
a supersymmetric generalization of the seesaw mechanism can accommodate massive 
neutrinos in SUSY models with conserved -R-parity. The superpotential of the Type-I 
seesaw extension of the MSSM is given by 


w = Wmssm + + -(Mu'UjNtm 


(19) 


where L* represents the chiral multiplet containing a SU{2)l lepton doublet (z/, 
and its corresponding superpartner, Hu represents the Y = 1 Higgs doublet and its 
Higgsino superpartner, Yf is a RH neutrino superheld, i,j are family indices, a,b are 
SU{2) indices and Cat is the antisymmetric SU{2) tensor.After EWSB, the neutrino 
mass matrix is given by ([^, as in the non-SUSY case. Moreover, the SUSY analogue of 
the Majorana mass term in the sneutrino sector leads to sneutrino-antisneutrino mixing, 
which could give rise to a same-sign dilepton signal at colliders 357- 361] . The neutrino 
Yukawa couplings lead to additional LFV effects in slepton masses, p —)■ ey in particular, 
through renormalization group effects in high-scale seesaw models 362 -368 . In low- 


scale SUSY seesaw models, new sources of LFV are present due to large neutrino Yukawa 


couplings and threshold effects from low-scale RH neutrinos and sneutrinos 369 372 


In such scenarios, the LFV rates of p — e conversion and /i —)■ 3e could be sizable, even 
bigger than the p —>■ ey rate. 

In low-scale SUSY seesaw models, the lightest superpartner of the RH neutrino 
with a small admixture of its LH counterpart could be another viable candidate for 
DM, if it happens to be the lightest supersymmetric particle (LSP) |354, 373 -377 


The sneutrino LSP scenario leads to distinct collider signatures, such as long missing 
transverse energy tail and enhanced same-sign dilepton signal in gluino and squark 
cascade decays 378-382], unlike the neutralino LSP in the MSSM scenario. 


6.2. R-Parity Violation 

Within the SM, the requirement of gauge invariance automatically guarantees B and 
L conservation for all renormalizable interactions. However, this is not the case in 
the general MSSM, and the following B and L violating terms are allowed in the 
superpotential: 


W, 


RPV 


= 5,6,,+ \,uea,L^L]El + + A'' 


jjL jNcrn r)cn 


( 20 ) 


where i,j,k are family indices, a, 6 are SU{2) indices, l,m,n are color indices, and 
the chiral multiplets respectively represent the {u,d)L,u‘i,d^j^,e‘i and the 

corresponding superpartners. Within the MSSM, these terms are forbidden by imposing 
an additional global symmetry that leads to the conservation of R-parity, given by 
R = (—1)3(s-D-2S S being the spin of the component held) 383 . The hrst term 
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in (20) leads to bilinear i?-parity violation (BRPV), while the remaining three terms 
collectively give rise to trilinear i?-parity violation (TRPV). 

The RPV models provide an alternative way to incorporate massive neutrinos with 
the minimal particle content of the MSSM (for a review, see e.g. 384| ). The BRPV 
model is the simplest one 44-^, and has several interesting consequences that can be 
probed at collider experiments (see e.g. 385- 387| ). Since the distinction between the 
matter doublet superhelds Lj and the Higgs doublet superhelds is lost in these 

models, it allows the mixing of neutrinos and neutralinos, sleptons and Higgs bosons, 
charged leptons and charginos. One linear combination of neutrino helds develops a 
Majorana mass at tree-level via mixing with Higgsinos, while the other combinations 
can acquire masses at loop-level via the trilinear couplings. Note that the trilinear 
coefficients A, A' and are constrained from data on various low-energy B- and L- 


violating processes 388 


The collider phenomenology of RPV models has quite distinct features from that 
In particular, the LSP is unstable, and hence, not all SUSY decay 


of the MSSM 384 


chains lead to a large missing energy at colliders. The phenomenology of pair-produced 
SUSY particles is also modihed due to new RPV decay chains. In addition, SUSY 
particles can now be singly produced, e.g. s-channel resonant production of sneutrinos 

and charged sleptons in hadron collisions 392f 394 . RPV 


m e 


~e collisions 389 391 


models also lead to sneutrino-antisneutrino mixing 395 and other low energy effects. 


such as 396-398 


7. Conclusions and Outlook 

The discovery of neutrino oscillations has provided us with the hrst conclusive 
experimental evidence for the existence of new physics beyond the SM. Therefore, 
just as the postulate of the very existence of the neutrino led to the formulation of 
the theory of weak interactions, an essential ingredient for the stupendous success of 
the SM, a clear understanding of the neutrino mass mechanism could as well be the 
hrst beacon of physics beyond the SM. Therefore, it is very important to explore the 
experimental signatures of various neutrino mass models to pin down the underlying 
new physics. In this brief review, we have discussed some low-scale neutrino mass 
mechanisms accessible to current and future experiments. In particular, we focused on 
the simplest Type-I seesaw model and summarized the current experimental constraints 
on the sterile neutrino mass and its mixing with active neutrinos. We have discussed 
the future discovery prospects of a heavy neutrino, within the minimal setup as well 
as involving extended gauge/Higgs sectors, with a particular emphasis on the energy 
frontier, in light of the upcoming run-H phase of the LHC and the proposed future 
colliders at both energy and intensity frontiers. A better picture of the neutrino portal 
might have far-reaching implications for the beyond SM scenarios in general, including 
the puzzles of matter-antimatter asymmetry and nature of dark matter in our Universe. 
In this context, we should emphasize the importance of complementary and synergetic 
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explorations in the low-energy sector at the intensity frontier, as well as cosmological 
observations at the cosmic frontier, a combination of which is essential to fully unravel 
the mysteries of the neutrino world. 
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